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In this paper we present a measure of quantum correlation for a multipartite system, defined 
as the sum of the correlations for all possible partitions. Our measure can be defined for quantum 
discord (QD), geometric quantum discord (GQD), or even by the entanglement of formation (EOF). 
For tripartite pure states, we show that the multipartite measure for the QD and the EOF are 
equivalent, which gives a way to compare the distribution and the robustness of these correlations 
in open quantum systems. We study dissipative dynamics for two distinct families of entanglement: 
a W state and a GHZ state. We show that, while for the W state the QD is more robust than the 
entanglement, for the GHZ state this is not true. It turns out that the initial genuine multipartite 
entanglement present in the GHZ state makes the EOF more robust than the QD. 
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I. INTRODUCTION 

Entanglement in composite quantum systems leads to 
many puzzling paradoxes in quantum theory The 
importance of entanglement is recognized unanimously 
but it is well known that even separable quantum states 
possess correlations that cannot be simulated by classi- 
cal systems; e.g., the nonclassical correlations captured 
by quantum discord (QD) For bipartite state.many 
efforts were done to detect and quantify QD 0, 0] and 
find the connection with entanglement [a |9|. Moreover, 
much attention is being devoted to the quantify of quan- 
tum correlation in multipartite systems [ToL M- Such 
measures help us to understand the distribution of quan- 
tum correlations and provide a way of studying dissipa- 
tive dynamics in many-part quantum systems. 

Although it is claimed that the QD is more robust 
than the entanglement in open quantum systems, little 
concrete evidence has been published to corroborate this 
claim where multipartite systems are concerned. Actu- 
ally, for two qubits there is a good deal of evidence that 
it is true |l2| , but can we extend this rule to multipartite 
quantum systems? For three qubits, for example, there 
exist two distinct families of entanglement fl3j : so, how 
robust, then, is the quantum discord for members of each 
family? Could entanglement be more robust than dis- 
cord, against noise arising from the environment, in mul- 
tipartite systems? There are very few studies on these 
questions in the literature 14| and that is the focus of 
this article. 

To investigate the robustness of the quantum corre- 
lations in open quantum systems for tripartite states, 



we define a multipartite measure of quantum correla- 
tion slightly different from that employed by Rulli and 
Sarandy [Hj . Those the authors define a measure of the 
global multipartite quantum discord as the maximum of 
the quantum correlations that exist among all possible 
bipartitions. Here, to attain an average measure, we de- 
fine global QD as the sum of correlations for all possi- 
ble bipartitions. Despite being only subtly different, our 
measure now accounts for how the quantum correlation 
is distributed in the tripartite system and certainly gives 
a better insight into the robustness of these correlations 
in open quantum systems. 

This paper is organized as follows: in section II we 
present the formal definition of the entanglement of for- 
mation, the usual quantum discord, and the geometric 
measure of quantum discord. In section III we define our 
multipartite measure of quantum correlations (MMQC) 
based on the sum of correlations for all possible bipar- 
titions. In section IV we present an analytical solution 
to the MMQC for a thrcc-qubit pure state and we show 
that for a general tripartite pure state the MMQC based 
on the usual quantum discord and on the entanglement 
of formation are equivalent. In section V we extend our 
analysis to a tripartite mixed state considering the dis- 
sipative dynamics of three qubits and in section VI we 
summarize our results. 



II. QUANTUM CORRELATIONS 

In this manuscript we consider three well-known mea- 
sures of quantum correlations: entanglement of forma- 
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tion, usual quantum discord, and the geometric measure 
of quantum discord. Below, we present the formal defi- 
nition of each one of these measurements. 



state with zero discord, i.e., 5^ B = 0, is a state whose 
information is not disturbed by local measurements; it is 
called as a classical-quantum (CQ) state. A CQ state is 
of the form [j| 



A. Entanglement of Formation 

The entanglement of formation is a measure of en- 
tanglement defined more than fifteen years ago by Ben- 
nett, DiVincenzo, Smolin, and Wootters [ll|. Although 
very different from QD, it is connected with the lat- 
ter a monogamic relation d, [l6| and has a nice opera- 
tional interpretation. It is defined as follows. Given a 
bipartite system A and B, consider all possible pure- 
state decompositions of the density matrix pab\ that 
is, all ensembles of states with probability pi such 

that pab — J2iPi\^i)(^i\- F° r eacri P ure state, the 
entanglement is defined as the von Neumann entropy 
of either of the two subsystems A and B, such that 
E(Hf) = Sa = S b , where S A ■= S(p A ) = -Tr(pA log/u), 
Pa being the partial trace over B, and there is an anal- 
ogous expression for Sb- The EOF for a mixed state is 
the average entanglement of the pure states minimized 
over all possible decompositions, i. e., 



E{ PAB ) = mm^piEi^i). 



(1) 



Whereas this is very hard to calculate for a general bipar- 
tite system, for two qubits there is an analytical solution 
given in the seminal William K. Wootters paper 



x = ]Tpfc|fc)(fc|®/°f 



(4) 
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where {pk} is a probability distribution, \k) an arbitrary 
orthonormal basis in subsystem A and p^ a set of arbi- 
trary density matrices in subsystem B. 

Denoting the set of all CQ states on Ha <8 H B as f2 0j 
then it is natural to think that, the further a state p 
is from f2o, the higher is its quantum discord. Indeed, 
we can use the distance from p to the nearest state in 
f2o as a measure of discord for the state, and this is the 
idea behind the geometric measure of quantum discord 
(GQD) [2l|. Thus, the GQD was introduced as: 



D{p) 
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where fJo denotes the set of zero-discord states and \ \X — 
Y\\ 2 — tr(X-Y) 2 stands for the squared Hilbert-Schmidt 
norm. Note that the maximum value reached by the 
GQD is | for two-qubit states, so it is appropriate to 
consider 2D as a measure of GQD hereafter, in order to 
compare it with other measures of correlation [22j . 

Interestingly, an explicit expression for GQD in a gen- 
eral two-qubit state can be written [2l|. In the Bloch 
representation, any two-qubit state p can be represented 
as follows: 



B. Usual Quantum Discord 

Quantum discord is a well known measure of quantum 
correlation defined by Ollivier and Zurek [j| about ten 
years ago. It is defined as 



' AB 



Iab — J 



AB 



(2) 



where Iab = Sa + Si 



and Jab ^ s the classical correlation [1 



Sab is the mutual information 
Explicitly, 



Jt~ R = max 
AB {n fc } 



sa - y^^PkS{p A \k) 



(3) 



where p A \k = Tr B (U k PAB^k) /Tt A b (Jl k pAB^k) is the 
locally post-measured state after obtaining the outcome 
k in B with probability pk ■ The quantum discord mea- 
sures the amount of the mutual information that is not 
accessible locally QilHiJ and generally is not symmetric, 
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XiO 1 ®I + y^ \y t l®a l 



i=l 



(6) 

where I is the identity matrix, a l (i = 1, 2, 3) are the three 
Pauli matrices, Xi = tr(cr I <8 T)p and yi = tr(I ® a l )p 
are the components of the local Bloch vectors x and y, 
respectively, and = tr(cr I (8cr- 7 )p are components of the 
correlation matrix T. Then the GQD of p is given by 



D{ P ) = \{\\x\ 



\T\ 



(7) 



A 



being the largest eigenvalue of the matrix K = 
xx* + TT* and ||T|| 2 = tr(TT*). The superscript t de- 
notes the transpose of vectors or matrices. Furthermore, 
it is important to mention that an analytical solution for 
a bipartite system of dimension 2 x N has been given in 

m. 



BA- 



C. Geometric Measure of Quantum Discord 



III. MULTIPARTITE MEASURE OF QUANTUM 
CORRELATIONS 



Intuitively, quantum discord can be viewed as a mea- 
sure of the minimum loss of correlation due to measure- 
ments in the sense of quantum mutual information. A 



Definition 1. For an arbitrary Af-partite state 
Pi-.-n-, the multipartite measure of quantum correlation 
Q (px-.-n) is defined as follows: 
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Let p be an iV-partite state, and p and v be any subsets 
among all possible partitions. The multipartite measure 
of quantum correlations (MMQC) is defined as the sum 
of the quantum correlations for all possible bipartitions, 

N 

Q(pi-n)= ( 8 ) 

where M^i^ is a measure of quantum correlation that 
can be given by the geometric quantum discord, the usual 
quantum discord or the entanglement of formation. Here, 
the subset between (•) is the measured one in the case of 
GQD or QD and can be ignored for the entanglement of 
formation. As we can see from the definition given in 
Eq. ([5]). our measure is symmetrical and, more impor- 
tantly, it gives zero if the state has just classical-classical 
correlations. 

To elucidate the MMQC defined above, let us consider 
the case of a tripartite state. Explicitly, 

2 (PABC) = M A(B) + M B (A) + Ma{ ) + M C (A) 

+ M B (C) + ^C(B) + Ma(BC) + M BC(A) 

+ M B ( AC ) + Mac{b) + Mc(ab) + Mab(c)- (9) 

IV. MMQC ANALYTICAL SOLUTION FOR 
THREE QUBITS PURE STATE 

Our starting point is an analytical solution to the 
MMQC for three-qubit pure states. For the case of ge- 
ometric measure of quantum discord, an analytical solu- 
tion can be found with the help of the result presented 
in [23|. Actually, from that result, an analytical solu- 
tion can be obtained for GQD even for three-qubit mixed 
states. The analytical solution for three qubits can also 
be obtained for the usual EOF. With the help of Wootters 
concurrence fl7j , each entanglement measurement involv- 
ing this sets of one subsystem (e. g. E B (c)) can be ob- 
tained trivially by calculating the eigenvalues of the den- 
sity matrix. Furthermore, we note that the entanglement 
measure involving a set of subsystem and a set of two (e.g. 
Ea(bc))i is given by the von Neumann entropy of one of 
the partitions. For example, Ea(bc) — Sa = S B c- 

Finally, to calculate the MMQC for the usual quantum 
discord, we note the result given in [25, where the au- 
thors show that the sum of the quantum discord for all 
possible bipartitions, involving sets of one subsystem, is 
equal to the sum of the entanglement of formation for all 
possible bipartitions, 

Ea(b) + E B (A) + Ea(c) + E C (A) + E B (c) + E C {B) = 

$A(B) + ^B{A) + &A(C) + ^C(A) + $B(C) + ^C(B)- (10) 

Moreover, since we are considering pure states, the quan- 
tum discord measure involving a set of one and a set of 
two subsystems is given by the entropy of one of the par- 
titions, exactly as the entanglement of formation. For 
example, E A (bc) = $a(bc) = S bc(A) = s a = S BC - 



Thus, for general tripartite pure states, the MMQC 
defined in Eq. (JSJ is identical for the QD and for the EOF 
which results in an analytical solution for the MMQC for 
three-qubit pure states, for the QD as well. To confirm 
this, we note the result given in Eq. (|10[) . which is valid 
irrespective of the system dimension of the subsystems. 
On the other hand, for general tripartite states (not three 
qubits), an analytical solution does not exist, either for 
the EOF or for the QD. 

V. MMQC FOR THREE-QUBIT MIXED 
STATES 

To calculate the MMQC for three-qubit mixed states, 
we limit ourselves to studying a rank-2 density matrix. 
In this case, as we show below, a simple strategy can 
be used. Since we are considering three qubits, to calcu- 
late the MMQC for sets of one subsystem (Ea(b), Ea(c)> 
Sa(b)> &A(cy eta ) * s trivial. In this case, the EOF can 
be calculated analytically by means of concurrence, and 
the quantum discord can be calculated numerically by us- 
ing positive-operator valued measurements (POVM) [24j . 
The question is: how can we calculate the MMQC for the 
terms that involves a set of one subsystem and a set of 
two, i.e. Ea(bc) ? E B{A c), $a(bc)-> S b(ac)> etc - ? Here > 
in contrast to the case of pure states, the von Neumann 
entropy does not give the right answer. The answer to 
this question is given by the monogamic relation between 
entanglement of formation and quantum discord 0, [l6j . 

Since we are considering a rank-2 density matrix, an 
extra subsystem E that purifies ABC is a two-level sys- 
tem. Thus, the idea is to calculate the quadripartite 
pure state pabce and use the monogamic relation. The 
monogamic relation implies that the EOF between two 
partitions is connected with the QD between one of the 
partition and the third one that purifies the pair. To elu- 
cidate this further, we present in detail the strategy used 
to calculate the QD and the EOF for the rank-2 tripartite 
mixed state pabc- Here, we show the way to calculate 
just the terms Ea(bc) &nd ^a(bc) > but the strategy is 
analogous for all the other terms involving a set of one 
subsystem and a set of two given in Eq. ©. For the 
EOF, the monogamic relation says that 

Ea(BC) = $A(E) + Sa\Ei (11) 

where A(E) involves two qubits, since ABC is of rank 2. 
In this case, we can calculate $a(e) numerically. For the 
quantum discord, on the other hand, the result can be 
reached by 

$A(BC) = E A(E) + S A \E, (12) 

where Ea(e) can be calculated by means of the concur- 
rence. Moreover, 5^ c ^ can be calculated numerically, 
since the measurements are just over one qubit, in this 
case, subsystem A. Thus, following the recipe above, it 
is straightforward to calculate the quantum correlations 
for all terms of Eq. © . 
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VI. MMQC IN OPEN QUANTUM SYSTEMS 

To calculate the MMQC in open quantum systems, 
we study two special situations: firstly, a three-qubit W- 
state subjected to independent amplitude-damping chan- 
nel and secondly, the GHZ state, with independent phase 
damping. The reason for this specific choice is that, for 
the whole of the dissipative process, we have a rank-2 
density matrix [25| . In this case, we can use the strategy 
explained above to calculate the MMQC for the usual 
QD and the EOF. 

The dissipative dynamics can be calculated straight- 
forwardly by means of Kraus operators [26] . Since we 
consider independent environments for each qubit, given 
an initial state for three qubits p(0), its evolution can be 
written as 

where the so-called Kraus operators E a ,p n = E a ® Ep® 
£ 7 satisfy £ Q:(3i7 E\^^E aS ^ = I for all t. The opera- 
tors E{ a y describe the one-qubit quantum channel effects. 
We first consider a W-state subjected to independent am- 
plitude damping. This damping describes the exchange 
of energy between the system and the environment and 
is described by the Kraus operators Eq = ^/p(a x +ia y )/2 
and Ei = diag(l, y/1 — p), where p = 1— e~ rt , T denoting 
the decay rate, and a x and a y are Pauli matrices. In Fig. 




time 

FIG. 1: (Color Online) the red curve (solid) shows the EOF 
dynamics, while the blue curve (dotted) shows QD and the 
cyan curve (traced) the GQD: initial W state subjected to 
independent amplitude-damping channels. 

((T|) we show the dissipative dynamics of the MMQC for 
an initial state given by \W) = (|100) + |010) + |001))/\/3. 
We see that in this situation the QD is actually more ro- 
bust than EOF and GQD. For a very short time, the EOF 
resists but it decays fast while QD maintains it robust- 
ness. This result corroborates the idea that the QD is 
more robust than the EOF in open quantum systems. In 



the second, and more important case, we consider an ini- 
tial GHZ state subjected to phase damping. The dephas- 
ing channel induces a loss of quantum coherence without 
any energy exchange. In this case the Kraus operators are 
given by Eq = diag(l, yjl — p) and Ei = diag(l, y/p). In 




time 



FIG. 2: (Color Online) the red curve (solid) shows the EOF 
dynamics, while the blue curve (dotted) shows QD, and the 
cyan curve (traced) represents the GQD. All of them for the 
initial GHZ state subjected to independent phase damping 
channels. 

Fig. © we show the dissipative dynamics of the MMQC 
for an initial state given by \GHZ) = (|000) + |lll))/ v / 2- 
Here, we see a very interesting result. Contrary to what 
is claimed in the literature, the EOF is in fact more ro- 
bust than the QD for this kind of initial condition and 
dissipative channel. EOF is sustained for a longer time 
than QD. 

This peculiar situation occurs because of the conserva- 
tive relation between the EOF and the QD Q. For the 
W state, the QD between the parts of the system is sus- 
tained and, consequently, these subsystems need to have 
more entanglement with the environment to preserve the 
conservative relation. Indeed, it is impossible to create 
or destroy some of the EOF (QD) without destroying or 
creating the same amount of QD (EOF). For the W state 
subjected to the amplitude-damping channel, each sub- 
system becomes more entangled with the environment 
but creates little discord. So, the QD is sustained in 
the system to maintain the conservative relation. On the 
other hand, for the GHZ state we have a genuine mul- 
tipartite quantum correlation, since the MMQC for the 
EOF or the QD is given by: 

Q (PABC) = Ma(BC) + M BC (A) + ^B(AG) 

+ M AC{B) + M C (ab) + M AB(C) . (14) 

To calculate Eq. (|14[) for the quantum discord, we use, 
as explained above, the monogamic relation between the 
EOF and the QD. In this case, each one of the six terms 
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is given by 



A(BC) 


— ^BC(A) 




\- Sa\e, 




'B(AC) 


= $AC{B) 








'C(AB) 


= ^AB(C) 


= E C (E) ~ 


\~ Sc\E- 


(15) 



However, there is no entanglement between each sub- 
system and the environment, i.e. £yi( e) = E B ^ = 
Ec(E) = 0, since the MMQC for the QD can be cal- 
culated analytically by means of the conditional entropy 
Indeed. 

Qqd(pabc) = 2(S a \e + S b \e + S c \e)- (16) 

For the EOF, the situation is a little different. In this 
case, each one of the six terms is given by 



Ea(bc) 


— E BC (A) 


— $A(E) - 


1- Sa\e, 




Eb(ac) 


= Eac(b) 


= ^B(E) ' 


1- Sb\e, 




Ec(AB) 


= E A B(C) 


= ^C(E) ~ 


\~ Sc\E- 


(17) 



The crucial difference between Eq. (fT5|) and Eq. (fl7|) 
is that, while each part of the system does not become 
entangled with the environment, it does create QD with 
it . Hence, S^ E ^ = <5^ B ) = &c(e) ^ an< ^' consequently, 

QEOF {PABC) = QQD (PABC) + 1{&A(E) + S B(E) + S C(E))- 

This is an important result and a direct consequence of 
the conservative relation between the EOF and the QD. 
For the GHZ state subjected to phase damping, there is 
no entanglement between each subsystem and the envi- 
ronment but there is QD. Thus, the QD that is created 
with the environment needs to be compensated by the 
entanglement retained in the system, making the EOF 
more robust than the QD in this particular situation. It 
must be emphasized that this is a direct consequence of 



the GHZ being a genuine multipartite entangled state, 
since this means that the EOF between sets of one party 
(Ea(b), Ea(e), E B (C), etc.) is always zero during the 
dissipative dynamics. 



VII. CONCLUSION 

In this paper, we have presented an alternative mea- 
sure of multipartite quantum correlations. Our measure 
gives a novel and intuitive means of comparing the ro- 
bustness of entanglement and discord in multipartite sys- 
tems, against the detrimental interaction with the en- 
vironment. We analyze two distinct initial conditions, 
which represent different families of multipartite entan- 
glement. We show that the robustness of the EOF and 
QD depends on the family of entanglement to which the 
initial state pertains, putting into question the superior 
robustness of QD in open quantum systems. Actually, 
for a three-qubit W state, QD proves to be more ro- 
bust, but the same can not be said of the GHZ state. 
We show that this behavior is related to the way the 
multipartite quantum state is quantum correlated with 
the environment. For the GHZ state subjected to inde- 
pendent phase-damping channels, contrary to the QD, 
the individual qubits do not get entangled with the en- 
vironment, and the initial multipartite entanglement is 
thus preserved for a longer time than the QD. We be- 
lieve that the discussion presented here may contribute 
further to the understanding of the distribution of entan- 
glement and discord in open quantum systems. 
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